Bell Function Values Approach to Topological Quantum Phase Transitions 
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We investigate the relation between Bell function values (BFV) of the reduced density matrix 
and the topological quantum phase transitions in the Kitaev-Castelnovo-Chamon model. We find 
that the first order derivative of BFV exhibits singular behavior at the critical point and we propose 
that it can serve as a good and convenient marker for the transition point. More interestingly, the 
value of the critical point can be analytically obtained in this approach. Since the BFV serves as a 
measure of nonlocality when it is greater than the classical bound of the correlation functions, our 
work has established a link between quantum nonlocality and phase transitions. 
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Topological phase of some strongly correlated quan- 
tum many body systems is a new kind of order that 
depends on the system topology It has attracted 
great interest recently because it can exhibit remarkable 
phenomena such as quasiparticles with anyonic statis- 
tics. An archetypal physical realization of such phase is 
in the quantum Hall system [pj, which bears many un- 
conventional characteristics, such as fractional statistical 
behaviors and ground state topological degeneracy that 
cannot be lifted by any local perturbations || . A partic- 
ular interest in topological ordered states is their robust- 
ness against local perturbations which can lead to several 
consequences such as topological insulators Q and topo- 
logical quantum computations || . 

Not surprisingly, the unconventional properties of 
topological phase might result in exotic critical phe- 
nomena, which cannot be characterized by the Landau- 
Ginzburg- Wilson spontaneous symmetry-breaking the- 
ory where the correlation function of local order param- 
eters plays a crucial role For example, the quan- 
tum phase transition between an Abelian and a non- 
Abelian topological phase in chiral spin liquid might be 
characterized by global flux and generalized topological 
entanglement entropy Q|. More remarkably, for time- 
reversal invariant anyonic quantum systems, Gils et at, 
have recently showed that the topological phases could be 
uniformly described in terms of fluctuations of the two- 
dimensional surfaces and their topological changes [^j. 
However, an universal characterization and detection of 
topological phase and its transitions still pose a big chal- 
lenge despite a vast amount of prominent works dealing 
with this problem. 

During the past few years, several important concepts 
in the quantum information field have been borrowed 
to characterize quantum phase transitions (QPTs) and 
topological quantum phase transitions (TQPTs), these 
including entanglement ||, fidelity fidelity suscep- 



tibility [[Ll| , and discord [[12|, etc. A brief review of the 
progress related to this issue is given in Ref. [|l3| and 
references therein. Notwithstanding the great successes 
in marking QPTs and TQPTs in some physical systems, 
each approach above has its own disadvantages Jl3[ . Take 
the fidelity approach for example, to witness the QPTs, 
one has to find out the exact ground state. However, 
for most of the physical systems, finding out the exact 
ground state is very difficult. In addition, it is also a 
challenge to measure the fidelity in experiment on scal- 
able systems. An alternative choice is to use Bell function 
values (BFV) as defined in expression (|^) below, which 
indicates the correlations of a quantum system and mea- 
sures quantum nonlocality when it is greater than the 
classical bound of the correlation functions Jl4| ]. Actu- 
ally, besides entanglement, quantum nonlcoality is also 
a central nonintuitivc phenomena of quantum mechan- 
ics and it plays a key role in many quantum information 
and computation processes, such as quantum key distri- 
bution (KQD) Jig] , nonlocal quantum computation [ [Ls| , 
etc. Naturely, one would ask whether nonlocality can 
mark QPTs and TQPTs? 

In this Letter, we propose the use of BFV as the marker 
of TQPTs and provide a positive answer to the above 
question by investigating the relation between BFV of 
the reduced density matrix and the TQPTs. The motiva- 
tion for choosing BFV is two-fold: (i) BFV can measure 
the nonlocality of a quantum system, thus it might estab- 
lish a connection between nonlocality and TQPTs, which 
belong to two different aspects; (ii) To get the BFV in 
an experiment scheme, one only has to do some measure- 
ment on the qubits instead of knowing exactly the ground 
state. Thus, our approach has its advantages in exper- 
imental schemes. The discussion here is mainly based 
on the Kitaev-Castelnovo-Chamon model |l7|, which ex- 
hibits a second-order TQPT at the critical point. Our 
results indicate that the first order derivative of BFV 
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shows singular behavior at the transition point. More in- 
terestingly, through this approach, one can analytically 
obtain the critical value of the transition point. Finally, 
using BFV to signal TQPTs and QPTs in other systems 
is also briefly discussed. 

Bell function values (BFV). — The famous Bcll- 
Clauser-Horne-Shimony-Holt (Bcll-CHSH) inequality for 
two entangled spin-1/2 (or qubit) particles, which has al- 
ways provided an excellent test-bed for experimental ver- 
ification of quantum mechanics against the predictions of 
local realism, is given by the inequality [JL8[ : 



1 = Qn + Q12 + Q21 - Q22 < 2, 



(1) 



where Q i6 = J r p(X)Xi (nf 1 , \)X 2 (n* 2 , \)d\ is the cor- 
relation function with Xk (n^ fc , A) denoting the m-th ob- 
servable on the fc-th particle (here i,j,k,m = 1,2); T 
is the total space of the hidden variable A and //(A) is 
a statistical distribution of A, satisfying J r p(X)d\ = 1. 
Quantum mechanically, the above inequality is violated 
by all pure entangled states of two qubits and the 
expression of the correlation function for any two-qubit 



state p reads: Q~< = Tr[(nf 1 • a) ® (n^ 2 • a)p\. Here n* fc 
(to, k = 1,2) are the unit vectors in three-dimensional 
Hilbcrt space and a is the Pauli matrix vector. Based on 
the Bell-CHSH inequality ([]]), the Bell function values is 
defined as: 



<M(p) 



(2) 



where I® = Qn + Q% + Q21 ~ Q22 an d the maximization 
is performed over all possible vectors n^ fc . Generally 
speaking, for every specific two-qubit quantum state p, 
we need to carry out the procedure of the maximization 
to obtain its BFV 98{p). Fortunately, in Ref. |(J, the 
authors introduced another method to calculate 33{p), 
which can circumvent the tedious maximization. It was 
proved there that 



&{p) = 2^ Vl + v 2 , (3) 

where v\ and v 2 are the two greater eigenvalues of the 
3x3 symmetric matrix Jzf^Jzfp; J*f p is a 3 x 3 matrix with 
elements defined by (-2 , p ) sr = Trf/oo^ ® o T \ (<j, r = 1, 2, 3) 
and Jfl is the transpose of Jzf p . In the experimental 
situation, in order to obtain the BFV, the observers of the 
first (second) qubit should carry out two measurements 
n^ 1 • ex and n 2 1 • a (n^' 2 • a and n^ 2 • a), just the same as 
in many Bell-CHSH inequality testing experiments [ pl| . 

The Kitaev-Castelnovo-Chamon (KCC) Model. — The 
physical model we consider in this article was introduced 
by Castelnovo and Chamon p7|, which is a deformation 



of the Kitaev toric code model |22 . The Hamiltonian of 
the KCC model with periodic boundary conditions reads: 

jr eT 2 VeT 2 y £T 2 
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FIG. 1: (Color online) An illustration of the Kitaev- 
Castelnovo-Chamon spin-lattice model and its map to the 
2D Ising model: In the KCC model, each black dot on the 
edge of the lattice represents a qubit and V and T denote 
the vertex and the face, respectively; The blue line across the 
lattices stands for a string operator along the non-trivial loop 
on the torus T 2 . The system is invariant under transforma- 
tion along the blue line |2~^ |. In its corresponding 2D Ising 
model, the qubits live on the vertex (green dots). Mapping: 
of = 6 r 6 r i, where i is the bond between the neighboring ver- 
tices (r, r'). Thus, for i and j nearest neighbors, the mapping 
gives (crfcr*) — {9 r 8 r i9 r i>8 r ) — (9 r iQ r ii), namely, the nearest 
qubits in KCC model become next-nearest in the correspond- 
ing 2D Ising model. 



where J m ^Je > 0, /3 is a coupling constant; A\> = 

Iljgv 17 j an d Bjr = Il : /ejF cr ! arc tnc vcr t cx an< i f acc 
operators in the original Kitaev toric code model p2| . A 
brief sketch of this model is shown in Fig. |l| The ground 
state in the topological sector containing the fully mag- 
netized state |0) = I \\ ■ ■ ■ f) can be analytically ob- 
tained ITmI: 



see 



i(ff)/2, 



10), 



(4) 



where Z{(3) = J2 g eg e ^ 1 ai ', Q is the Abelian group 
generated by the vertex operators {Ay} and o-j(g) is the 
value of spin at site j in state g\0). Obviously, when 
f3 = 0, \G{fi)) reduces to the topologically ordered ground 
state of the toric code model [ p2| , while when f3 — > 00, 
\G(j3)) becomes the fully magnetized state |0). At the 
point j3 c = i ln(V2+l) there exists a second-order TQPT 
where the topological entanglement entropy Stopo = 1 for 
(3 < f3 c changes to S top o = for j3 > (3 C @ . 

As shown by Castelnovo and Chamon, there exists a 
one-to-two mapping between the configurations {g} = Q 
the configurations {9} of the classical 2D Ising model [Il7| . 
In the mapping, the Hamiltonian of the Ising model has 
the form Rising = — c ^J2(rr') ®r®r' , where ^ is a cou- 
pling constant and 6 r ,6 r > = +1 or —1 depending on 
whether or not the corresponding vertex operator Ay 
is acting on the site r. Thus af = 6 r 6 r : with i being 
the edge between the nearest neighboring vertexes. An 
illustration of this mapping is shown in Fig. [j]. 

Signaling TQPTs by BFV. — Since the BFV introduced 
in expression (|^) only account for two-qubit states, we 
need to calculate the reduced density matrix of two qubit 
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FIG. 2: (Color online) Numerical results of the first order 
derivative of BFV £i§(pij) versus p for different 5j3. (a) i and 
j axe the nearest case; (b) i and j are the next-to-nearest case. 



Pij based on the ground state |G(/3)) and the symmetry 
of the Hamiltonian Jif. It was shown in Ref. |23 that pa 
has the following form (details are given in Ref. |2J] 
and references there in): 



Pij = -[I + (of )« + o-|) + K0-|Kcr|], 



(5) 



where I is the 4x4 identity matrix. Based on the 
Eq. (§), the BFV can be calculated by using the sim- 
plified formula for 38{pij) in Eq. (||). For conve- 
nience and simplicity, we only concentrate on two cases 
where i and j are nearest and next-to-ncarcst neigh- 
bors, respectively. In the thermodynamic limit, the 
mapping to the 2D Ising model gives that (erf) = 
(Mr') = -coth(20)[7r + (4tanh 2 (2/3) - 2) (x)]/(2tt), 



where JT^) = ^ dp(l 
2sinh(2/?)/cosh 2 (2/3). For the calculation of (erfcrj), the 
equivalence between the 2D Ising model and the quantum 
ID XY model yields: 

(i) For i and j the nearest case, (of er|) = (8 r 8 r >) = 
i J£ d(j){[smh~ 2 (2/3) — cos 0] cos — sin 2 (jn}/{sm 2 4> + 
[sinh~ 2 (2/3) — cos 0] 2 } 1 / 2 . Summarizing all the relations 
above enables us to obtain the BFV BS(pij). For the i 
and j nearest case, the numerical results for the first order 
derivative of BFV d "^^' J ^ are displayed in Fig. [§(a), from 



X 2 sin 2 tp) : / 2 and \ 



which we see a distinct rapid increase of 



dp 



round 



the critical point /3 C ps 0.44. Note that we only focus on a 
small neighboring region of (3 around the TQPT point /3 C , 
namely 0.4 < (3 < 0.5. The smaller the 8/3, the greater 
the rapid increase. When 5/3 — > 0, d "^' J ^ — >• +oo, indi- 
cating its singularity at the TQPT point /3 C . 

More interestingly, after long tedious but straightfor- 
ward calculations, we arrive at an analytical formula for 
d ^dp ' J ' wn i cn can enable us to obtain the analytical 
value of /3 C : 



~dj3 



ksch 2 (2/3) sin 2 0T(0, 0)/v, 



(6) 



Here T(0,/3) 8 coth (2/3)csch 2 (2/3)/(l - 

2 cos Resell 2 (2/3) + csch 4 (2/3)) 3 / 2 . What is interesting is 
that we can analytically obtain the critical point from 



the Eq. (g). To this end, one can rewrite T((f>,(3) as 
T(0,/3) = 2\/2coth(2/3)/{[i(csch 2 (2/3) + l/csch 2 (2^))- 
cosc/>] 3 / 2 csch(2/3)}. Obviously, Eq. (0) has only one sin- 
gular point because csch 2 (2/3) + l/csch 2 (2/3) > 2 and so 
the singularity happens at csch 2 (2/3) + l/csch 2 (2/3) = 2, 
namely (3 C = iln(\/2 + 1). This explicitly exhibits one 
of the advantages of BFV approach to TQPTs. 

For i and j the next-to-nearest case, direct calcula- 
tions show that (ofo-J) = cosh 2 {(3* ) (T 2 1 - T_ 2 T ) - 
sinh 2 (/?*)(7? - T 2 r ), where T K = ± ft d0[(£ - 
cos (/>) cos K<p + 7 sin <f> sin K^]/[(7sin 0) 2 + (£— cos 0) 2 ] 1 / 2 . 
Here tanh(/3*) = e~ 2 ^ 7 = [cosh^/?*)]" 1 and £ = 
(1 — 7 2 ) 1 / 2 / tanh(2/3). We also plot the numerical re- 



sults of 



dPJ(pi 
dp 



111 Fig. |(b). From this fi gure, one can 



observe that the rf ^^' J ^ has a singularity at the TQPT 

point (3 C . It is also obvious that dSg ^ 3 - behaves quite 
similarly between the nearest and next-to-nearest case. 
This result accords with the results in Ref. ^3|, where 
the reduced fidelity and reduced fidelity susceptibilities 
are only slightly different between the two cases, respec- 
tively. Since BFV indicates the qubit correlations in the 
system, it seems that the correlations of nearest qubits 
is similar as that of next-to-nearest qubits in topologi- 
cal ordered states. This is, to some extent, counterin- 
tuitive because in many physical systems the interaction 
between nearest particles is usually greater than that of 
next nearest particles. 

It is worthwhile to note that the reduced density ma- 
trix p^ is diagonal. That indicates the correlations be- 
tween any two local spins in the ground state of KCC 
model are always classical. Consequently, the Bell func- 
tion values of pij cannot be greater than 2, the classical 
bound. Another interesting consideration here is similar 
to Ref. pi| : we can calculate the BFV between a lo- 
cal qubit denoted by i and the rest of the whole lattice 
by rewriting the ground state as \G(fi)) = W+\S?)\§)i + 
where ^ ± 2 = (1 ± (0 o>o 0o,i))A \&) and \2) 
are two orthogonal normalized vectors. Consequently, we 
can regard |G(/3)) as a simple pure two-qubit entangled 
state. In this case, the BFV has a one-to-one monotonous 
relation with entanglement |2^] , thus also with quantum 
discord since for pure two-qubit state, the quantum dis- 
cord is the same as entanglement of entropy [ p6| . As 
a result, the BFV should behave similarly as the quan- 
tum discord does at the critical point (3 C (for details, see 
Ref. p3j). However, it is worthwhile to clarify that there 
is a distinctive difference between the BFV approach and 
the quantum discord approach. For the quantum discord, 
its value becomes trivially for the reduced two-qubit 
state, thus cannot signal the TQPT at the critical point. 
Nevertheless, as shown in the former paragraphs, the first 
order derivative of BFV is an excellent marker of the 
transitions. In addition, the physical meanings of BFV 
and quantum discord are different. Generally speaking, 



quantum discord is a measurement of the 'quantumness ' 
of a system. While, BFV measures the nonlocality of the 
system when it is greater than the classical bound. In this 

case, @(\G(P))) 
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2 i yi + 4^ + 2 ^_ 2 > 2. Thus the BFV 

3§(\G((5))) can measure the nonlocality of ground state. 
This establishes a new link between quantum nonlocality 
and TQPTs. 

Summary and Discussion. — To summarize, based on 
the KCC model, which exhibits a second-order TQPT 
at the critical point, we have introduced the BFV ap- 
proach to TQPTs. Our results show that BFV serves as 
an accurate marker of the transitions. Since the BFV 
also serves as a measure of nonlocality, which is a pure 
quantum phenomenon and cannot be described by any 
local realism theory, our work has established a new 
link between quantum nonlocality and phase transitions. 
Furthermore, experimentally, this approach only involves 
two measurements on two qubit, thus it might be more 
convenient to implement in experimental schemes. Actu- 
ally, the optical lattices and trapped ions might provide 
suitable experimental test-bed for our results j27). 

What is also notable is that this approach is applica- 
ble to other models. For instance, for the model recently 
introduced by Son et al., which is described by a clus- 
ter Hamiltonian Jt(X) = - Eili(°"f-i°"f of+i+^M+i) 
and has an exotic phase transition at the critical point 
A = 1 ^8|, our numerical results show that the first or- 
der derivative of BFV can explicitly capture the transi- 
tion. The cluster Hamiltonian above can be simulated 
in a triangular configuration of an optical lattice of two 
atomic species |28), thus also leading to the possibility 
of testing the BFV approach experimentally. To investi- 
gate the BFV in QPTs, we have also considered the one 
dimensional Ising model and XY model. The numerical 
results show that the first order derivative of BFV ex- 
hibits singular behavior at the critical point, too. Thus 
this approach is useful for both QPTs and TQPTs. 

It would be interesting and significant to apply this 
approach to QPTs and TQPTs in various physical sys- 
tems, such as quantum spin Hall system, both theoret- 
ically and experimentally. It would also be interesting 
to use BFV based on other Bell inequalities to inves- 
tigate QPTs and TQPTs. More specifically, studying 
the BFV of the pure ground states based on multipartite 
Bell inequalities, such as the famous Mcrmin-Ardchali- 
Belinskii-Klyshko (MABK) inequality might shed 
light on the behavior of the quantum nonlocality of the 
whole system in QPTs and TQPTs. 
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